STABILITY FOR THE MULTI-DIMENSIONAL BORG-LEVINSON THEOREM 

WITH PARTIAL SPECTRAL DATA 
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Abstract. We prove a stability estimate related to the multi-dimensional Borg-Levinson theorem of de- 
termining a potential from spectral data: the Dirichlet eigenvalues Aj. and the normal derivatives d4>k/dv 
of the eigenfunctions on the boundary of a bounded domain. The estimate is of Holder type, and we allow 
finitely many eigenvalues and normal derivatives to be unknown. We also show that if the spectral data is 
known asymptotically only, up to 0{k~'^) with a 2> 1, then we still have Holder stability. 

1. Introduction 

In 1988, Nachman, Sylvester, Uhlmann [17] proved an rt-dimensional version of the classical one dimen- 
sional Borg [B] and Levinson |14) theorem: one can determine uniquely a potential q in the Schrodinger 
operator —A -I- g, from knowledge of the Dirichlet eigenvalues and the traces of the normal derivatives of 
the normalized eigenfunctions on the boundary of a bounded domain. The proof is based on relating the 
spectral data to the Dirichlet-to-Ncumann map Ag_;^ for all frequencies A, see (|3.4p below. Then we recover 
the potential by applying Aq_A to high-frequency solutions. In fact, one frequency is enough for uniqueness 
P5] but the recovery then is only logarithmically stable [T1[TH]. The link between the spectral data and Ag_A 
was also noticed by Novikov in [T5] . 

The first stability estimate for this problem of conditional Holder type was proved by Alessandrini and 
Sylvester [2] . A variant of this result was given by the first author in [7j and extended later by Bellassoued and 
Dos Santos Ferreira [5] to the case of the Schrodinger operator on a simple Riemannian compact manifold with 
boundary. The main idea in the approach initiated by Alessandrini and Sylvester consists in transforming 
the stability estimate for the spectral problem into a stability estimate for the inverse problem of determining 
the potential in a wave equation from the corresponding hyperbolic Dirichlet-to- Neumann map Hq. They 
gave an explicit formula relating h-q-\ to Hq. On the other hand, reconstruction of a potential from Hq can 
be done in a stable way using geometric optics f5UJ [52] . A powerful method based on the boundary control 
method was initiated by Belishev and developed by Belishev, Katchalov, Kurylev, Lassas and others. We 
refer to [TT] for more details. This method proves uniqueness of recovery of the coefficients of a general 
elliptic operator from the corresponding hyperbolic Dirichlet-to-Neumann map H . An important ingredient 
of this method is the unique continuation result of Tataru [24l [25] which makes it unlikely to provide Holder 
stability when the latter might hold. In the metric case, conditions on the metric are needed apparently, see 
e.g., [2I1[5]. 

In [3], a logarithmic type stability estimate with a partial hyperbolic Dirichlet-to-Neumann map Hq was 
proved, under the assumption that the potential is known near the boundary. The proof of this result relies 
on a qualitative estimate of continuation from boundary data for the wave equation. The result for the wave 
equation yields a log-log type stability estimate when the traces of the normal derivatives of the normalized 
eigenfunctions are known only on a part of the boundary. Recently, an extension of this result to a log type 
stability estimate was proved in 0] under an additional assumption in terms of the X-ray transform of the 
potential. The general problem of determining the potential from Dirichlet eigenvalues and the traces of the 
normal derivatives of the normalized eigenfunctions on a part the boundary is still an open problem. An 
earlier work by Isakov and Sun [8] proves stability estimates for a partial hyperbolic Dirichlet-to-Neumann 
map in the special case where the measurements are made in the intersection of the boundary of the domain 
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with a half space. They estabhsh a Holder type stability estimate in dimension three and a logarithmic type 
stability estimate in dimension two. 

An inverse spectral problem with different spectral data was considered by Kurylev, Lassas and Weder 
[12]. The case of a singular potential was considered by Paivarinta and Serov [19], following the approach 
by Nachman, Sylvester and Uhlmann [17] . 

In [9] 110] , Isozaki proved that if we are missing a finite number of eigenvalues and eigenfunctions traces, 
then we can still determine q uniquely. He also mentioned that in fact, a sharp enough asymptotic formula 
would be enough, see also Theorem 12.11 below. 

The purpose of this paper is to prove that even if a finite number of spectral data is missing, we still 
have conditional Holder stability, see Theorem 12.21 In fact, in Theorem 12.31 we prove something more: if 
the spectral data are known only asymptotically, with a sharp enough estimate of the remainder, similarly 
to Theorem 12. 1[ then the asymptotic data determine the potential in a Holder stable way. Theorem 12.31 
combines the previous two theorems but it has stronger assumptions. 

Our work is inspired by the method introduced by Isozaki in [9] . It is built on high frequency asymptotics 
techniques. The main advantage of this approach is that it is a direct method. In other words, it is not 
necessary to relate the spectral problem to an hyperbolic Dirichlet-to-Neumann map for the wave equation. 
Moreover, the results we obtain this way are stronger than the preceding ones. 

2. Main Results 

Let be a bounded domain of M" that we assume, for simplicity, of class C°° . Its boundary will be 
denoted by F. li q G L°°(r2), we denote by A{q) the unbounded operator acting on L^(51) as follows 

A{q) = -A + q, DiA) = H^{n)nH^in). 

We recall that the spectrum of A{q) consists in a sequence of eigenvalues, counted according to their multi- 
plicity. This sequence can be ordered in the following way : 

-oo < \i{q) < X2{q) < ...< Xk{q) < • . . ^ +oo. 

In the sequel we will use the notation X{q) = {Xk{q)}. We note that as an immediate consequence of the 
classical min-max principle, we have X{q) G A(0) + £°°, where i°° is the usual Banach space of bounded 
sequences equipped with its natural norm. We denote by {(f>k{Q)} an orthonormal set of eigenfunctions, each 
one related to Afc(q). Note that (f>ki<l) is defined only up to a factor of modulus 1 when Xk{q) is a simple 
eigenvalue; and when Afc(g) is multiple, then we have more freedom, and the natural way is to think about 
eigenspaces. Within each such eigenspace there are infinitely many choices of eigenfunctions and to get the 
best results in the theorems involving their traces on F, in (12. 1[) . we should minimize (take infimum) over all 
such possible choices. 

Our first result is a uniqueness theorem under the assumption that the spectral data are asymptotically 
"very close" . As a partial case, we recover the result in [9j about uniqueness with a finite number of spectral 
data missing. It was noted in [5^ that such a result should hold. 

Theorem 2.1 (Uniqueness). Let qi,2 G L°°{n). Let, for some A > 0, and all k — 1,2 ... , 

\Xk{qi)-Xk{q2)\<Ak~", a>l, 

(2.1) 1 

WdMiqi) - a.0fc(<Z2))||L^(r) < Ak~^, ^ > ^ " 

Then qi = (72 • 

We shall use the following useful upper and lower bounds for eigenvalues. Let M > be given. Then 
there exist two c* > and c* > 0, depending only on M and fl such that, for all q £ L°°(il) satisfying 

l|g|U°=(n) < M, 

(2.2) c*fc2/" < Afe(g) < c*fc2/", fc > 1. 

From the usual elliptic regularity estimate we have : for any e > 0, there exists a constant Cg 

(2.3) \\d.Mq)\\L-(r) < \\d.MQ)\\HHr) < CAk{qf^^^'^^ 
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for any q G L°°{il) satisfying |1<z||l«>(o) < M. Here the constant C depends only on fl and M. 
Let 

(2.4) m>n/2 + 3/4 

be a fixed integer. Using estimate (|2.2|) in ()2.3|) . we easily obtain that the sequence {k~^"^^^^\\di,ipk{q)\\ L^(r)} 
belongs to the Banach space of sequences such that the corresponding series are absolutely convergent. 
Since L'^{T) is a Banach space, this is equivalent to say that the series ^ k~'^^"-^^dij(pk{q) converges in L^{T). 
Let > be a fixed integer, we set X{q) — {Xk+N{q)} and 

So{qi,q2) = WHqi) - Hq2)\\i^- 

We consider the quantity 

Si{qi,q2) = ^kr'^"'^"'\\d^(pk+N{qi) - 9^</?fe+Ar(g2)||L2(r). 
fc>i 

Our second main result is the following stability theorem. 

Theorem 2.2 (Stability). Let q\, q2 G L°°{^) such that qi — q-z E i?o(f^) and 

lki||L~(n) + ll92lU-(n) + < M. 

Then there exists C = C{n, Q, m, M) > and < 7 = "f{n) < 1 such that 

hi - q2\\mn} < cs'^, 

where S = (5o(9i,92) + (5i (91,92) • 

We will see in the end of the proof of this theorem that we have an explicit value of 7 as a function of n. 

Next theorem combines the results of the previous two in one, see also Remark 12.11 It shows that for 
stable recovery, we only need to know that the data are asymptotically close. In particular, missing a finite 
number would not affect the stability but it will affect the constants, of course. 

Theorem 2.3. Let gi, (72 £ L°°(il) such that q qi — q2 ^ H^iyt) and 

IklllL-(O) + ||g2||L-(0) + \\q\\Hl(iT.) < M. 
Fix m satisfying ^2.4\ l- Let, for some 5 > Q, A > Q, 

(2.5) \\k{qi) - Mq2)\ < S + Ak~", k-^-^/^'+^Wd^Mqi) ~ duMq2)\\L^r) < s + Ak-^, 

with a > (4m — l)/(2n). Then there exists C — C{n, fl, m. A, a, M) > and < 7 = j{n, a) < 1 such that 

\\qi - '72||l2(o) < CS'^. 

Remark 2.1. Theorem 12.31 implies somewhat weaker versions of Theorem 12. II and Theorem 12.21 with a and 
/? in Theorem 12.11 required to satisfy stronger assumptions; and with a stronger norm of the traces of the 
normal derivatives of the eigenfunctions in Theorem 12.21 

Acknowledgments: We thank Gunther Uhlmann for his suggestions during the preparation of this 
paper. 

3. Preliminaries 

First, we consider a family of Dirichlet to Neumann maps parametrized by the spectral parameter. Let 
q G L°"{il) and A G p{A{q)), the resolvent set of A{q). Following well known results on existence and a priori 
estimate (e.g. |15j), for each / G H^/^{T), the following boundary value problem 

( {-A + q- X)u = inn 
^ ' ' \ u = f on F 

has a unique solution u{q,X){f) G H^{V,) and / — > di,u{q, \){f) defines a bounded operator from i7^/^(F) 
into H^^/^(r). We denote this operator by A{q, A). It extends to a mcromorphic family with poles at the 
eigenvalues. 
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Let 92 e L^{n), A e p{A{q,)) n p(Afe)) and / e i/i/2(r). Then u = ufe, A)(/) - wfe, A)(/) is the 
solution of the foUowing boundary value problem 



(3.2) 



{-A + qi - X)u = {q2 - qi)u{q2, X){f ) in 
u = on r 



Therefore, according to the classical a priori estimate, A{qi,X) — A{q2,X) defines a bounded operator 
from i?i/2(r) into H^/^{r). 

Next, A{q,X), q £ L°°{n), is symmetric as form on C°°{T) x C°°(r) with respect to the duality form 
{f,9) Jr fgdSa:. Indeed, 

{A{q,X)f,g)^ [ (GAF-FAG) dx ^ [ {{q - X) - {q - X))FG dx = 0, 
Jn Jn 

where F = u{q,X){f) and G = u{q,X){g). Since A{qi,X) — A(g2, A) is a bounded operator on H^^'^{T), its 
transpose, which is the same operator, is bounded on H^^/^(r). By interpolation, 

(3.3) A(gi,A)-A(g2,A) :i/^(r)^i/^(r), |s| < 1/2 

is bounded, as well. Note that for q smooth, A(q,X) is a pseudo-differential operator of order 1, see |13) . 
while one can see that the difference p.3p is of order —1 either by (|3.2p or by calculating the first few terms 
of the symbol in the spirit of I13|. Then, in particular, (13. 3p can be improved. 

The following formal representation of A{q, A) providing a relationship between the spectral data and the 
family of DN maps A{q, A), appears in [17] 

(3.4) A{q,X)f J2^—^^d,^^{q){f,d.^,{q))mr). 

The series on the right hand side is not absolutely convergent in some special cases, at least, even if considered 
as a form. A possible way to justify it, suggested in [2], is to show that some high order formal derivative 
converges. Set 

A('"n'Z,A) :=_aM(5,a). 

It was then shown in [171 [7] that for m ^ 1, the series converges absolutely because Afc(q)^™^^ ~ ^-2(m+i)/ra 
decays rapidly when m ^ 1, while the traces of the eigenfunction on F have a fixed polynomial bound. 

Lemma 3.1. Let q e L°°{i}), (p{q) an orthonormal basis, f G H^/'^(r), m > n/2 + 3/4 and X G p{A{q)). 
Then 

(3.5) A('")(g, A)/ = -ml ^ ^^^^^^^^ ^yn+i ( / fdu^l) d.Ml), 

where the series converges absolutely in iJ^/^(r) and therefore in £^(r)0 

We can adapt Lemma 2.32 in jTj which goes back to [T7] to the complex case and with H^^^iT) in place 
of ij3/2(r). We obtain 

Lemma 3.2. Under the assumptions of Theorem \2.'A for any positive integer I and < e < 1/2, there exists 
a constant > 0, that can depend only on M , f2, I and e, such that 

(3.6) ||A(-'")(gi,A)-A«(g2,A)||^(^v.(r).L^(r))<^^^, 0<j<Z, 5RA < -2Af, = 



^ In view of the proof of Lemme 2.28 in [7], we can prove, using a density argument, that the result in Lemma 13. 1 1 remains 
valid for / e L'^{T). 
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The importance of this lemma is to provide some decay of the difference A{qi, A) — A((72, A), as SRA — > — oo, 
in suitable norms; and similarly for the derivatives. Such an estimate is not surprising, of course, because 
when 5RA — — cx), A moves away from the spectrum. 

By Lemma \3l2\ 



(3.7) A(gi,A)-Afe,A) 



dAi 



dXr 



(A('")(gi,A„0-A(™)(g2,A„)) , A 



in C{H ' {T),L (r)), where, for non-real A, the integrals above are taken over the lines 9Aj = 9A, j = 
1, . . . , m — 1. The estimate p.6p in Lemma [3.21 for j > 1, shows that each integral is absolutely convergent; 
and the same estimate for j > shows that the initial condition after each integration is zero at — oo. Now, 
plugging (|3.5p into p.7|) provides a direct formula expressing the difference of the DN maps in terms of the 
spectral data. On the other hand, integrating term by term is not justified, which is the main reason to work 
with the differentiated series (13.51) . 

At the end of this section, we recall a lemma in j9] which is the basis for the proofs of the main theorems. 

Let (p\^uj{x) = e*^"'^, A G C \ (— oo, 0], with the standard choice of the branch of the square root, and 
u) € S""-^, we consider 



S{q){X,UJ,9) 



A(<7,A)((^A,^)^A,_(,da(x), A e p{A{q)) \ (-oo,0], uj, 0e§" 



Following Lemma 2.2 in f9| we have, for A G p{A{q)) \ (-oo,0], uj, 9 e §" 



(3.8) 



Siq,X,io,e) = --\9-uj\ 



dx 



+ 



e-*^('-")-"g(x)da: - / Riq,X){qipx,^)qipx,.edx, 



where R{q, A) — {A{q) — A) ^ is the resolvent. 

We fix ^ G S"-i and r] G §"-\ rj^C For t > 1, let 

9r = CrT] + 
UJr = Crf] - 

\/A7 = t + i. 

Let qi, q2 satisfy the assumptions of Theorem l2.2l and let (p{qi), (p{q2) be an arbitrary orthonormal basis. 
We fix < e < 1/2 and we set a = = (1 - 2e)/4. 

In the sequel C is a generic constant that can depend only on n, ft, M and e. Also, for simplicity, we 
drop the subscript in Ar, oJt and 9r- 

Using the classical estimate for the resolvent, where q — qi or ^2, 

\\R{q,z)\\cimn)) < z<^R, 



we easily prove 
(3.9) 



R{qi, X){qiipxM)qi^x,^edx 



R{q2,X){q2^px,Lo)q2V\-edx 



< 



C 



We deduce from identity (|3.8p and estimate p.9p . where the extension by zero outside w of g is still 
denoted by q, 



(3.10) 



l(9l-92)(e + -OI <- + \S{quX,9,Lo)~S{q2,X,0,Lo)\. 

T T 



Since qi — q2 is an entire function, uniqueness would follow if we can show that the difference in the right 
hand side tends to 0, as t — > oo. For stability, we need to estimate the same difference in terms of the 
spectral data. 
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4. Proof of the uniqueness result 
In this section, we prove Theorem 12. II 

Proof of Theorem \2.1\ Notice first that 

(4.1) 0A,^(a;) =e(*^-i)"-" =^ |0a,c.(x)| < C, xe^. 

By (1231), the difference A(qi, A) - A (92, A) is bounded in L^(r). 
We get from (jSlOl) . 

(4.2) |(gi - <72)((1 + */r)OI < C||A(gi, (r + zf) - A(q2, (r + *)')||£(L^(r)) + C/r. 
where C depends only on an a priori bound for HgiHiooi-f^-) + ||'Z2||l°°(si) and Vt. As a corollary, if 

(4.3) ||A(qi,(T + z)2)-A(q2,(T + i)')ll£(L2(r)) =0(1), asr->oo, 
then qi =52- 

Let us see what assumptions would imply (j4.3p . Let / and g be fixed with 

ll/llL2(r) = 1, ll.9llL2(r) = 1- 

We have, formally, 

((A(gi, A) - A(g2, A))/, g) = ^ (<zi), g)(a,<^,(<Zi), /) 



(4.4) 



where 



E T3T-^(^'>'^j(92),ff)(5.0,(g2),/) = /i + I2 + h, 



^ A-Aj(g2) 



/./ ^ E - 3^^) (5.0.(.r),..)(5.0,(<z.), /), 



3 

^2 



From now on, A = (t + i) . For /i, we get 
Note that for any e > 0, 

\{d.^,J)\ < ||9.0,||L^(r)ll/IU^(r) < CAj^^/^+^y^. 
We used the following estimates here: 

(4.5) \\d,(bj\\mr) < 115.0, ||^./2(r) < C||0,|1^3me/. < CeAf < c'^^^+'^f^. 

Therefore, 
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For A fixed and not among the eigenvalues, the series converges when —a + (3 + 2e) /n~A/n< —1, i.e., when 
a > \ — 1/n (one can always find e ^ 1 so that —a + (3 + 2e)/n — 4/n < —1). For A = (r + i)^, we apply 
Lemma 14^1 below, see also Remark 14.21 with ^ = ~a + (3 + 2e)/n, < e ^ 1 and v = 2 to conclude that 

|/i| < Cr(i-")"+2^ 

when a < 1 + 1/n + 2e/n. Choose a — 1 + 3e/n, e < 1/n to get 0{t~'^) above. 
To estimate I2, write 

■(3/2+e)/n 

(4.7) |/2|<CE IX X r ^1 • 

^ |A-Aj(g2)| 

We apply Lemma |4T1 again, with v = \ and /i = — /?+ (3/2 + e)/n, see also Remark 14.31 Since /3 > 1 — l/(2n), 
we can always find < e ^ 1 so that /3 > 1 — l/(2n) + 4e/n. With that choice of /3, we apply the second or 
the third inequality in (|4.9p . Note that < 2/n — 1; in fact, ^ < 2/n — 1 — 3e/n, and {fj, + l)n/2 < 1 — 3e/2, 
see ()4.9p . If, in addition, — 1 < /i, we apply the second inequality to get 

I/2I < C'^r'+(^+^)"/2"'' < Cr"'/^ 

If ^ < — 1, we apply the third inequality in (14. 9p that gives us the even better estimate I2 = 0{t~^). 
We treat I3 in a similar way. 

Therefore, (|4.3I) is satisfied, so (71(C) = 92(0 for any ^. Therefore, = q2. □ 

Remark. Let us discuss briefly why formula (4.4) is valid under assumption (2.1). To this end we set 
/(A) = Ii + 12 + Hi where /i, I2 and /a are as above. Then one can prove in a straightforward way that / 
is analytic in p{A{qi)) H p{A[q2)). By using Lemma 3.1 and the fact that weak analyticity is equivalent to 
strong analyticity, we deduce, where m>n/2 + lis fixed, 

/(™)(A) = ((A(™)(<Zi,A)-A(")(g2,A))/,5), A e p(A(gi)) n ^(^(92)). 

Therefore 

rn—l 

/(A) = ((A(<zi, A) - A(<Z2, A))/, g) + Y, afeA^ A e p{A{q^)) n p{A{q2)). 

k=l 

From Lemma 3.2 we know that 

hm ((A(qi,(T + z)2)-A(g2,(T + *)2))/,5) =0. 

On the other hand, from the proof of Theorem 2.1, we easily see that 

lim I{{T + if) = 0. 

r— >+oo 

Hence at = 0, < fc < m — 1. 

The following lemma was used in the proof. Its proof, on the other hand is based on Lemma 14.21 below. 

Lemma 4.1. Let 

00 

Then for p, < 2z//rt — 1 and X = t + i, t > 0, the series converges absolutely and for any e > 0, t > 1, 

Ct^^+i)"-!-"^ for 2/71 - 1 < p, 

(4.9) |/((t + i)2)| < { CeT'+{f^+l)n/2-^ j^^ -1 < p < 2 / 71 - 1, 

Ct-"^ for p<-\ . 

Remark 4.1. The condition 2/n — \ < p is not compatible with the convergence condition p < 2v jn — 1 
when ^ < V < \. In that case, only the second and the third inequalities above are applicable. The case 
V = \ was of particular importance above, see ()4.7p and also Remark 14.31 
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Proof. We will sketch first the idea. Since Xj{q) ~ i^^", we get 
(4.10) 



„2/n 



4t2)' 



The series and the integral are convergent, if — 2i'/n < —1. Make the change x = t to get 

roo ,b J. 

(4.11) I' = C 6:=(M + l)n/2-l. 

Jl {{t ~ T^Y + At^) ' 



The convergence condition takes the form h — v < —1. By Lemma 14.21 below, for any e > 0, and for 
-1 + 2v/n, 

(4.12) p ^(j^(i^+i)n-i-v for < (a* + l)n/2 - 1, i.e., for 2/n - 1 < ^, 

(4.13) /» < c'^^E+(M+iW2-.^ for -1 < (^ + l)n/2- 1< 0, i.e., for -1 < /i < 2/n - 1, 

(4.14) /" < C,T~'' for {jl + l)n/2 - 1< -1, i.e., for /i < -1 . 

We proceed to the actual proof now. The only step that needs to by justified is (|4.10p above. Let N{r) be 

1/2 

the counting function of the square roots A^- [q] of the eigenvalues \j{q) counted with their multiplicities, 

1/2 

i.e., N{r) = #{Xj {q) < r}. We will use the Weyl asymptotic formula with a sharp estimate of the reminder 
term (reference?) 

(4.15) A^(r) =c„r" + 0(r"-i). 
Another way to write this asymptotic formula is 

(4.16) Ay'(g)=c~„ji/« + 0(i). 
As a consequence, 

(4.17) |Ay'(9) <A„ 

with some A„ > 0. Moreover, An can be chosen independent of q as long as q belongs to a fixed ball in 
We start with 

~ jfj, <^ -fj. 



((r2 - 1 - A,(<z))2 + Ar^'^ U ((^' ^ " ^M)^ + - 1)) 



Set Ti = Vt^ — 1, r > 1. Write for simplicity / :— rf + 1 + i)^) to get, after replacing ri by r again, 



1^1 <E 



rt((r2-A,(g))2+4r2))''/2- 



Clearly, it is enough to prove the estimates (|4.12P " (|4.14p for /. Split the sum above into three parts; the 
first, /i, over the terms with A^^^(g) < r — 2A„; and the second, I2, for t — 2 An < A^^^(g) < t + 2^„; and 
the third containing all terms with Xy^{q) > t + 2 An- When Xy^{q) < r — 2A„, we have 

(4.18) |t2 - X^(q)\ = r2 - X,{q) > ~ (c„ji/" + Anf. 

On the other hand, Cnj^^" + An < Ay^(g) + 2 An < t, therefore, the term on the right above is non-negative. 
Then 

(4.19) (r^ - A,(g))2 > (r^ - (5„//" + AJ^)' . 
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When X/ (q) > r + 2A„, we have 

(4.20) |r2 - A,(g)| = \,{q) - > (£„//" - A.^f - t\ 

Similarly, Cni^^"^ — > ^^^(9) ^ 2A„ > r, therefore, the term on the right above is non- negative and we 
can take squares of both sides of the inequality to get 

(4.21) (r2 - A,(g))2 > (r^ - (5„ji/" - A^f)\ 
Therefore, 



h < 



Since the r.h.s. of (|4.18p is non-negative, each summand the first series in (I4.22p is an increasing function of 
j. Therefore, it can be estimated by above by the integral 

A r v!:dy 

Ji ((r2-(5„yi/« + A„)2)2 + 4T2)-/2- 

The change c„y^/" + An = x^^" reduces the analysis of /{ to that of (|4.11[) . Then we use Lemma [4.21 to 
estimate /f, to conclude that the estimates (j4.9p in Lemma HTTl applv to /J, and therefore to /i, as well. 

We analyze /a in a similar way to conclude that the estimates (j4.9p apply to /a, as well. 

The analysis of I2 is straightforward. In the interval [r — 2A„,t + 2A„] there are 0(r"^-'^) square roots 
of eigenvalues Aj's. The numerator admits the estimate j'' < CA"^^^(q) < C'r"''. The denominator can be 
estimated by below by Ct^ . Therefore, 

h < Ct"-i+"^-'^, 

which is as in (j4.12p . and stronger than (I4.13p . (|4.14p . This completes the proof of the lemma. □ 

Remark 4.2. The only property of the eigenvalues that we used was the estimate (|4.17l) . The lemma remains 
true if we replace (A — Xj{q)y by the product (A — Xj{qi)Y^{X — Xj{q2)y^, vi + V2 = where v\,2 > 0, 
ly = + 1^2- 

We used the following lemma in the proof. 

Lemma 4.2. Let P be as in ^JJ^. Then Let b — v < —1. Then, for all e > 0, 

(4.23) /« < Cr^^'+i-'^, for < b, 

(4.24) /« < Ct'^+'^+I"'', for-l<b<0, 

(4.25) I^<C,T-'', forb<-l. 

Remark 4.3. Computer algebra calculations in a few special cases show that (|4.23p is sharp at least for 1^ — 2, 
while ()4.24p . (|4.25p are not when 1/ = 1; the estimate in case (|4.24p then seems to hold as in (|4.23p . with a 
possible logarithmic term. The only case in the proof of the theorem above when b is not strictly positive 
is when v — \ and 6 < 0, small; and we want to have — o(l). We cannot apply (14.231) then, and we use 
(I4.24p . The loss of the factor in (|4.23l) compared to (|4.24l) for < -6 < 1 then is not essential. 

Proof. Perform another change of variables t = ts -\- r'^ . Assume in what follows that t > 1. We get 

^ ' ~ Jr-^-r [r's^ + 4t2)-/2 " (,,2 + 4)'^/2 

For s > r^^ — r, and 6 > 0, we have 

(4.27) (./r + l)''<(|.s|/r+l)''<(|.| + l)^ 
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When b < 0, we have s/t + 1 > t ^, therefore, 

(4.28) (s/r + l)^<T2|''l. 
Therefore, by g^Z]), 

(4.29) /» < T^'+'-' I + ds < Cr^b+i-. for < 6 < - 1. 
On the other hand, for 6 < we get by (|4.28p . 

(4.30) < Ct^-" for & < 0, 1< v. 

This is weaker than what we need to prove (|4.24p . To get the stronger estimates in the lemma, we 
apply the Holder inequality to (|nT|) as follows. Write t'' = t^'H''^ with &i,2 < 0, b = bi + 62- Then, for 
1/pi + l/p2 = 1, Pi > 1, P2 > 1, we have 



The first integral is the same as but with b — pibi, D — piiy. Since b — v < —1, and pi > 1, we have 
Pib — piv < —1, therefore that integral converges again. By (|4.30p . if, £> > 1, it is 0(t^^''). If p2|fe2| > 1, 
then the second integral in (|4.3ip is convergent. Therefore, 

(4.32) /» < Cr^i^Pi^^/Pi = Ct^/p'^". 

It remains to see what choice of pi , p2 , bi, 62 satisfying the conditions above gives the best decay. We need 
to minimize (up to an e > error) the factor 1/pi. Choose 61 — eb, 62 = (1 ^ ^)b, with < e ^ 1. We need 
to satisfy the conditions 

P2I62I > 1, i>^pii^>l, Pi,2 > 1, — + — = 1. 

Pi P2 

Set l/p2 = (l-2e)|6|. If e < 1/2, and \b\ < 1, we havep2 > 1 as required. Then 1/pi = l-(l-2e)|6| > 0. 
We also have pi > 1. Then 1/pi — = 1 — (1 — 2e)|6| — v can take any value greater than 1— |6| — z/ = 1 + b — iy. 
The latter is negative by assumption, so we can take 1/pi — < when e <C 1; then v > I. This proves 
(|4.24p . The proof of (|4.25p is immediate since (in view of the definition of /") the denominator is bounded 
by below by (4t^)^/^, while the numerator is integrable over (1, 00). □ 

In next section, we would need an estimate of I{z), as 5Rz < 0. 

Lemma 4.3. Let /(A) be as in ^J^. Then for /i < 2i'/n — 1, v > and 3A < the series converges 
absolutely and for 3 A < — 1, 

(4.33) |/(A)| < c\^x\^^'+^^"/^-^-'' < C. 
Proof For QA < 0, we have |A ~ Xj\ > |5A| + Xj. Therefore, for SA < -1, 

' ^ ^' (|3A| + A,r ^ (|3A|+c,j2/«)'' Jo (|3A|+c.a;2/")^- 

Make the change of variables |3A|t = .t^/" to get easily 



|/(A)| < C|3A|(^+i)"/2-i-", 2vln < -1. 



□ 



Remark 14.21 applies to this case as well. 
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5. Proof of the stability estimate 



We start with ((XTU| . Set q = qi - q2- We have 



T 1 7 c 

m)\<m+-o\ + - sup |vg(^ + -e)-ei. 

T T 0<S<1 T 



On the other hand, 
with some c = c(f2). Hence, 

(5.1) m)\<m+k)\ + ^e'^\\qh^^n)- 

Until the rest of the proof, we choose X — t + i, t > 1. Combining (|3.10p and (|5.1I) yields 

c|g(OI<- + 

T 

Integrate this in the ball |f | < r" to get 



Cm)\ < i + ^e'^ + |5(gi,A,0,L.)-5(<Z2,A,0,w)|. 

T T 



2 - T 



c \mvdi < 



.an ^a(2+^^) 



\i\<r- 



t2 r2 



+r""|5(gi,A,(?,c.)-5(g2,A,0,c.)|2. 

Choose a = 1/(2 + n) to get 

(5.2) c[ |g(Oprfe< -+r"/('^+')|5(gi,A,0,c.)-5(<72,A,^,c.)|2. 

J|5|<Tl/{2+r.) T 

Since g G iJi(M"), we have 

Iklli2(n) = \\q\\l2(^n) 

< / m)\'dc+-^jf-:^ f mmi'd^ 



< 



|^|<^i/(,>+2) ■ ■ r2/("+2) j|^|>^i/(„+2) 

\^iO\^d^ + —1— -llgll^i 

|^|<^l/(,^ + 2) r"^/l." + ^j ^ ' 



Then it follows from ()5.2p 

(5.3) C||g||i.(,,) < -5^1^ + r"/("+2) |5(qi. A, 0, lo) ~ Siq^, A, 

We now estimate the last term in the above inequality in terms of ||A(qi, A) — A{q2, A)||. We have 
\S{qi,X, 0, ^) - S{q2, X,9,Lu)\^\J^ <^A,-e[A(gi, A) - A(g2, A)]^a.c. da 
< C{n)\\A{qi,X) - A(q2, A)|| ||'^A,c^||//i/2(r)||'^A -0||L2(r). 
Here, and in the rest of the proof, || . || denotes the norm in C{H^/^(r), L^{T)). Since 

||</5A,a;||_H-i/2(r) < Cr^/^, \\(px-e\\L^(r) < C, 
we deduce easily from the last inequality 

|5(gi. A, 9, u) - S{q2, A, e,w)\< Ct1/2 II A(gi, A) - A(q2, A)|| . 

By (IS31), 

(5.4) qi<zlli2(o) < ^7^+r"/("+^)+^||A(gi,A)- A(g2,A)|p. 
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For q — qi or q2, we decompose A{q, A) in the following form 

A(g,A)=A(g,A)+A(g,A), 

where, for / e H^/^iT), 

(5.5) A(g,A)/ = ^ \ (]-x { I fduMl)do{x))^k{q)) e R-^I^T) (formally), 

(5.6) A(g, A)/ - > ' , , : , ( / fd.^k {q)da{x) ) d.^k {q) 



,^A,(g)-AV7p 



If c = c(il, M) is a constant such that Afc(g) < ck'^^^, for all fc > 1, then in a straightforward manner we 
obtain 

(5.7) l|A(-'")(9,^)ll < > 2ciV2/„ < 



In particular, ()5.7p implies 

(5.8) llA(g,A)|l < r>To, 

for some tq > 1 depending on and fl. 
In view of (|5.4p . estimate (|5.8p implies 

^lklli^(o) < 757^^:5) + ,2(n+3!/(„+2) + A(gi, A) - A(,„ A)|r 

Hence, 

(5.9) C\\q\\h^^^ < -^j^ + r"/(«+2)+i|| A(gi, A) - ^q^, A)|p. 

Let m — [71/2 + 3/4] + 1, where [71/2 + 3/4] is the entire part of [n/2 + 3/4]. We will use the estimate 
from Lemma 13.21 More precisely, we have 



(5.10) ||AW(gi,z)-A(^)(g2,z)|| < A G C, 5Rz < -2Af, < j < m. 

Here cr = (1 - 2e)/4 for some fixed e, < e < 1/2. 

It follows from (|5.7p that this estimate holds true if we replace A{qi, ■) by A{qi, ■). That is we have 

(5.11) ||A(^)(gi,z)-A«(g2,;2)|| < ^ ^ C, 3?A < -2Af, < j < m. 

In the sequel, changing tq if necessary, we can assume that Tq — 1 > 2M. For p > 29fiA, we set A = — p + A. 
From Taylor's formula, we have, for q = qi or (72, 

(5.12) A(,, A) ^ AW (,, A) + /' (i^in^_^A W (,, A + ,s(A - A))d.. 

^ A:! 7o (m-1)! 

We introduce the following temporary notations 

n?,A)='x:'(^^A(^)(g,A) 

r (^-;)"'(-^ -^)'" A(")(g,A + .(A-A))d.. 
Jo (rn-l)! 

Since SRA < — 2Af, a straightforward application of inequality (15. lip leads to 

(5.13) ||p(,„A)-7'(g2,A)||<^<-^. 
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By Lemma |3. 11 we know that 

In the sequel /i = /i(s) = A + s(A — A) = A — (1 — s)p and 

N{\) = min{/c > N: \k+i{q) > 23?A}. 
For 3?A » 1, we decompose A(™)(g,^)/ as foUows 

A(") {q, Ai)/ = a1™) (g, fi)f + A*") (q, m)/, 



where 



a1™^('7,m)/= XI 7m I fc>u^k{q)da{x))d,,ipk{q), 



We have 
where 



A;>A'(A) ^ '"''^'^ ^'^ 

a1™^(9iM)/ - A^™^ (52, A^)/ = /i + /2 + /a, 

JV(A) ^ 1 1 / /■ 

^ ,^5.1 Vxk{qr)-,r^^ - (A.(,2)-M)"-^Ky/^''^'=('^)'^-(-))^''^'=('^^)' 
"^^^^ 1 / /• X 



k=N+l 
JV(A) 



-?^3/= E TT ^—^-z^i [ fdu'Pk{q2)dcr{x))\d^(pk{qi) - d^(pk{q2) 

For /3 > 4/n + 1 fixed, we get easily 



2 

k=N+l 



As ||¥'fc(gi)lli2(r) < Cfc*/", we have 



Therefore 



Similarly, we obtain 



Setting 
we deduce 



JV(A) 

Y: fc-^|!^.(gi)|!i.(r)<E^"'^'^"- 
fc=Ar+i fc>i 

\\h\\<C^^Jo{qi,q2). 

ll^2il + l|/3H<C y^,,,^, ^l(gi,g2). 
6 = (5o((7i,(72) +^1(91,92), 



Then the choice of /3 = 2(m + 5/4)/n, satisfying /3 > 4/ri + 1 by our choice of to, yields 
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(5.14) ||A("')fa,^)/-A("')fe,M)|| <C ^'^'jj^^^,^^^ S. 
Since 

\Xk{q) - Ml > Afc(q) - 5RA + (1 - s)p > Afe(g) - 5RA > 
we can proceed similarly to the case SRA < —2M (see the proof of Lemma 2.33 in |7], p. 72). We find 

l|A^^(9i,A')/-A^'"'((Z2,A^)/|| <C5. 

This and (|5.14p imply 

llA(™Hgi,M)/- A('")(g2,A^)/|| < C7V(A)2("+5/4)/"5. 
From the definition of A^(A), we have 

CiV(A)^/" < 23?A. 

Hence, 

IIA^^H'Zi^A*)/- A(™H(Z2,/i)|| < C3?A™+5/4^, 

and then 

l|A(™)(gi, /i)/ - A('")(<72, < Ct2(™+5/4)j. 

Therefore 

(5.15) ||7^(gl, A) - 7e(<z2, A)|| < Cp^T^f^+s/^)^. 
It follows from ^J^, ([ET^ and IjbAd^ that 

C7||A(qi,/i) - ~A{q2,fi)\\ < -1 +p"V2("+5/4)^. 

Plug this estimate in (15.91) to get 

(5.16) < ;^57^ +r^("+^^/("+^^(^ +P^"r^("+^/^)^ 

Let us choose p = with k = l/(2cr)Q 

This choice of p in ()5.16p yields 

(5.17) C^lklli^(n) < ^37^ + r2("+i)/("+2)+4™+4(™+V4)^2^ 
A standard minimization argument, with respect to r, leads 

C\\q\\mn)<S\ 

with 

1 



rt + 2 + 2(n + 2)(Km + m + 5/4) ' 

The proof is then complete. 



■^Note that since t > 1, p > 25iA is satisfied. 
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6. Proof of Theorem 12.31 
The starting point is (lOl) . Choose N = N{S) so that iV"" = 6. Then 

(6.1) \Xj{qi)~Xj{q2)\<6 + Ar" ^N-^ + Aj-'^ <{1+A)j-^ for j < TV. 
On the other hand, 

(6.2) |A,(gi) - A,(<Z2)| < S + Aj-" < (1 + A)5 for j > N. 

We have analogous inequahties for the expression in (j2.5p involving the eigenfunctions because the right-hand 
sides of the two inequahties in (|2.5p are the same. More precisely, 

(6.3) J-^"'^"+'\\d.Mli) - 5.0,('?2)||L^(r) < (1 + A)j-" for j < TV, 
i.e., 

(6.4) I15.0,(gi) - 9.<^,(g2)||L^(r) < (1 + A)r"+'"/"-' for j < N. 
Also, 

(6.5) r II a.<^,(<Zi) - 9.0,(g2)|U.(r) < (1 + A)<5 for j > iV. 

Now, define A(A) and A(A) as in (|53)) . ((CT with that A^. While A^ = A^((5) depends on (5, the upper 
bounds in (16. ip - (|6.5p do not. 

Estimates (|6.1I) and (|6.4I) are of the type (|2.1I) with the same a and with fi — a ~ 2m jn + 1. By the 
assumption on a, we have /3 > 1 — l/(2n), as required in (|2.ip . Next, the assumption on a together with 
(|2.4p imply a > 1 + 1/n; therefore, the first condition in (|2.ip holds as well. Therefore, by Theorem 12. H and 
its proof, see (|4.3p . 

(6.6) ||A(<zi,(r + z)2)-A(<Z2,(r + z)2)|| <Cr-^ r > 1, 

with some > depending on a. The constant C is independent of (5 because the number of terms in the 
sum (|5.6p does not matter — we can just complete it to an infinite series by adding zero terms — but the 
upper bound (|6.ip is independent of 8. Notice that in (|4.3p . the operator norm is in £(L^(r), L^(r)), which 
is stronger than the L{H'^I'^{^),L?{V)) norm which we use in this proof, denoted simply by || • ||. With 
A = r + i, (|6.6p replaces estimate (|5.8p but we note that in (j6.6p . we have the difference of two A's, instead 
of estimating each one. 

Proceeding as in the proof of Theorem 12. 2[ we combine (|5.4p and (|6.6p to get 

(6.7) q|g|li.(n) < ;^+r"/("+2)+i||A(gi,A)-A(g2,A)|p, 6, := min(0, 2/(n + 2)), 
compare to (|5.9p . 

To get the estimate (|5.1ip . we need an estimate replacing (15. 7p first. We will derive an estimate similar to 
(|5.7p for 9z < (only) but for the difference of two A's, as above. Indeed, following the proof of Theorem l2.11 
let us estimate first /i, see (gH) for 3A < 0. By LemmaSJl h < C\'^X\-'^+^ < C|9A|-i by choosing e in 
(|4.6p small enough. This proves that the contribution of Ii to (|5.7p is 0(|5A|~^), when j = 0, as claimed. 
When j > 1, we apply the same argument to the differentiated series, similarly to the proof of Theorem l2.2l 
Each differentiation increases ly in the application of Lemma 14. 3 1 bv 1, and we get 0(|5A|~^~^ ), as stated. To 
estimate I2, see (|4.7p . we apply Lemma again with /i = —/3 + (3/2 + e)/n (and /3 — a), and v = 1. The 
requirement (/i + l)n/2 < v is satisfied for e <C 1, and then we get I2 < C|SA|"^. As before, differentiating 
A(A), and applying the same argument, we see that each derivative brings another power of |9A|~^. To 
summarize, we have 

(6.8) ||A«(gi,^)-A«(g2,z)||<^^^, 3iz<Q. 
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Now, using (|6.8|) (instead of (15. yi) ). we prove (15.111) in the current setup. The rest of the proof is the same 

N<k<N(\) 



with J2k=^+i considered as J2N<k<N{\) because N may not be an integer anymore. In (|5.16|) . we get 



(6.9) < ^ + ^2(n+l)/(«+2)+4K™+4(m+2)^2^ 

The proof foUows by a minimization, as before. 
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